Statistical mechanics of 2D turbulence with a prior 
vorticity distribution 



P.H. Chavanis 

Laboratoire de Physique Theorique 
Universite Paul Sabatier 
118, route de Narbonne 
31062 Toulouse, France 



Abstract 

We adapt the formalism of the statistical theory of 2D turbulence in the case where the 
Casimir constraints are replaced by the specification of a prior vorticity distribution. A new 
relaxation equation is obtained for the evolution of the coarse-grained vorticity. It can be used 
as a thermodynamical parametrization of forced 2D turbulence (determined by the prior) , or 
as a numerical algorithm to construct arbitrary nonlinearly dynamically stable stationary 
solutions of the 2D Euler equation. 

Two-dimensional incompressible flows with high Reynolds numbers are described by the 
2D Euler equations 

^+U'Vw = 0, u = -Aip, u = -z x Vtb, (1) 
at 

where uj is the vorticity and ip the streamf unction. The 2D Euler equations are known to 
develop a complicated mixing process which ultimately leads to the emergence of a large-scale 
coherent structure, typically a jet or a vortex [1]. Jovian atmosphere shows a wide diversity of 
structures: Jupiter's great red spot, white ovals, brown barges,... One question of fundamental 
interest is to understand and predict the structure and the stability of these equilibrium states. 
To that purpose, Miller [2] and Robert &: Sommeria [3] have proposed a statistical mechanics 
of the 2D Euler equation. The idea is to replace the deterministic description of the flow 
uj(r, t) by a probabilistic description where p(r,a,t) gives the density probability of finding 
the vorticity level uj = a in r at time t. The observed (coarse-grained) vorticity field is then 
expressed as uj(r,t) = J pada. To apply the statistical theory, one must first specify the 
constraints attached to the 2D Euler equation. The circulation r = J uJdr and the energy 
E = \ J uJipdr will be called robust constraints because they can be expressed in terms of 
the coarse-grained field uj (the energy of the fluctuations can be neglected). These integrals 
can be calculated at any time from the coarse-grained field uj(t, t) and they are conserved 
by the dynamics. By contrast, the Casimir invariants // = / f(uj)dr, or equivalently the 
fine-grained moments of the vorticity = J uj n dr = J pa n dadr, will be called fragile 

constraints because they must be expressed in terms of the fine-grained vorticity. Indeed, the 
moments of the coarse-grained vorticity T^'^ = Juj n dr are not conserved since uj n ^ uf 1 (part 
of the coarse-grained moments goes into fine-grained fluctuations). Therefore, the moments 
must be calculated from the fine-grained field uj(t, t) or from the initial conditions, i.e. 
before the vorticity has mixed. Since we often do not know the initial conditions nor the 
fine-grained field, the Casimir invariants often appear as "hidden constraints" [4]. 
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The statistical theory of Miller-Robert-Sommeria is based on two assumptions: (i) it is 
assumed that we know the initial conditions (or equivalently the value of all the Casimirs) 
in detail (ii) it is assumed that mixing is efficient and that the evolution is ergodic so that 
the system will reach at equilibrium the most probable (most mixed) state. Within these 
assumptions 1 , the statistical equilibrium state of the 2D Euler equation is obtained by maxi- 
mizing the mixing entropy 

S[p} = - J p\np drda, (2) 

at fixed energy E and circulation V (robust constraints) and fixed fine-grained moments 
^n>i (fragile constraints). This optimization principle is solved by introducing Lagrange 
multipliers, writing the first order variations as 

5S - I35E -a5r-J2 a n ST f n 9 - = 0. (3) 

n>l 

In the approach of Miller-Robert-Sommeria, it is assumed that the system is strictly 
described by the 2D Euler equation so that the conservation of all the Casimirs has to be 
taken into account. However, in geophysical situations, the flows are forced and dissipated 
at small scales (due to convection in the jovian atmosphere) so that the conservation of the 
Casimirs is destroyed. Ellis et al. [7] have proposed to treat these situations by fixing the 
conjugate variables a n> \ instead of the fragile moments If we view the vorticity levels 

as species of particles, this is similar to fixing the chemical potentials instead of the total 
number of particles in each species. Therefore, the idea is to treat the fragile constraints 
canonically, whereas the robust constraints are still treated microcanonically. This point of 
view has been further developed in Chavanis [8]. The relevant thermodynamical potential 
is obtained from the mixing entropy (2) by using a Legendre transform with respect to the 
fragile constraints [8]: 

S x = S-^a n T^. (4) 

n>l 

Expliciting the fine-grained moments, we obtain the relative entropy 

drda, (5) 
where we have defined the prior vorticity distribution 

X (a) = exp{-5>„<A (6) 

^ n>l > 

We shall assume that this function is imposed by the small-scale forcing. Assuming ergodicity, 
the statistical equilibrium state is now obtained by maximizing the relative entropy S x at fixed 
energy E and circulation V (no other constraints). The conservation of the Casimirs has been 

1 Some attempts have been proposed to go beyond the assumptions of the statistical theory. For example, 
Chavanis & Sommeria [5] consider a strong mixing limit in which only the first moments of the vorticity are 
relevant instead of the whole set of Casimirs. On the other hand, Chavanis & Sommeria [6] introduce the 
concept of maximum entropy bubbles (or restricted equilibrium states) in order to account for situations where 
the evolution of the flow is not ergodic in the whole available domain but only in a subdomain. 
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replaced by the specification of the prior xO 7 )- Writing SS X — (35E — aST = 0, and accounting 
for the normalization condition / pdo = 1, we get the Gibbs state 

1 /*+oo 

p(r,a) = — x (a)e-^ +a > with Z= X (a)e^ +a > da. (7) 

^H r J J-oc 

This is the product of a universal Boltzmann factor by a non-universal function X (a) fixed 
by the forcing. The coarse-grained vorticity is given by 

57 = = Fm + a) with F ^=-^y^ w 

where x($) = x^e"^^. It is easy to show that = -u; 2 ($) < 0, where 

W2 = — cU 2 > is the local centered variance of the vorticity, so that F is a decreasing 
function [9]. Therefore, the statistical theory predicts that the coarse-grained vorticity uJ = 
f(ip) is a stationary solution of the 2D Euler equation and that the uJ — ip relationship is 
a monotonic function which is increasing at negative temperatures (5 < and decreasing 
at positive temperatures (3 > since Lo'(ip) = —j3u>2- We also note that the most probable 
vorticity (c)(r) of the distribution (7) is given by [10]: 

(a) = [(lnx) / r 1 (/3V' + a), (9) 

provided (lnx)"((c)) < 0. This is also a stationary solution of the 2D Euler equation which 
usually differs from the average value a;(r) of the distribution (7) except when X (a) is gaussian. 
We note that the oJ— relationship predicted by the statistical theory can take a wide diversity 
of forms (non-Boltzmannian) depending on the prior x( c7 )- The coarse-grained vorticity (8) 
can be viewed as a sort of superstatistics as it is expressed as a superposition of Boltzmann 
factors (on the fine-grained scale) weighted by a non-universal function x(c) [4]. Furthermore, 
the coarse-grained vorticity (8) maximizes a generalized entropy (in tJ-space) of the form [11]: 

Sp] = - J C(uJ)dr, (10) 

at fixed circulation and energy (robust constraints). Writing 5S — (35 E — a5T = leading 
to C'(lJ) = —Pip — a and uj'(tp) = —/3/C"(uj), and comparing with Eq. (8), we find that C 
is a convex function (C" > 0) determined by the prior x(°") encoding the small-scale forcing 
according to the relation [4]: 

C(uJ) = -j U F~\x)dx = - j" [{hi X )']-\-x)dx. (11) 

The preceding relations are also valid in the approach of Miller-Robert-Sommeria except that 
x(<t) is determined a posteriori from the initial conditions by relating the Lagrange multipliers 
a n> \ to the Casimir constraints r{'^ 1 . In this case of freely evolving flows, the generalized 
entropy (10) depends on the initial conditions, while in the case of forced flows considered 
here, it is intrinsically fixed by the prior vorticity distribution. 

In that context, it is possible to propose a thermodynamical parameterization of 2D 
forced turbulence in the form of a relaxation equation that conserves circulation and energy 
(robust constraints) and that increases the generalized entropy (10) fixed by the prior x( cr )- 
This equation can be obtained from a generalized Maximum Entropy Production (MEP) 
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principle in iU-space [11] by writing the coarse-grained 2D Euler equation in the form DuJ/Dt = 
— V • um = — V • J and determining the optimal current J which maximizes the rate of entropy 
production S = — J C"(lu)J • VuJdr at fixed energy E = J J ■ Vipdr = 0, assuming that 
the energy of fluctuations J 2 /2uJ is bounded. According to this principle, we find that the 
coarse-grained vorticity evolves according to [11, 8]: 

^ + u-VuJ = V-|d[vuJ+^^V^}, uJ=-A</>, (12) 

/ DVuj ■ Vj>d 2 r 1/2 1 

«*>— JD g^ r • = 7^W (13) 

where j3{i) is a Lagrange multiplier enforcing the energy constraint E = at any time. These 
equations increase the entropy (^/-theorem S > 0) provided that D > 0, until the equilibrium 
state (8) is reached. The diffusion coefficient D is not determined by the MEP but it can be 
obtained from a Taylor's type argument leading to expression (13)-b [8]. This diffusion coeffi- 
cient, related to the strength of the fluctuations, can "freeze" the relaxation in a sub-region of 
space ("bubble") and account for incomplete relaxation and lack of ergodicity [12, 13]. The re- 
laxation equation (12) belongs to the class of generalized Fokker-Planck equations introduced 
in Chavanis [11]. This relaxation equation conserves only the robust constraints (circulation 
and energy) and increases the generalized entropy (11) fixed by the prior vorticity distribution 
x(o~)- It differs from the relaxation equations proposed by Robert & Sommeria [14] for freely 
evolving flows which conserve all the constraints of the 2D Euler equation (including all the 
Casimirs) and increase the mixing entropy (2). In Eqs. (12)-(13), the specification of the 
prior x(c) (determined by the small-scale forcing) replaces the specification of the Casimirs 
(determined by the initial conditions). However, in both models, the robust constraints E 
and r are treated microcanonically (i.e. they are rigorously conserved). Furthermore, in the 
two- levels case oj G {a"o,<7i}, the two approaches are formally equivalent and they amount to 
maximizing a generalized entropy (10) similar to the Fermi-Dirac entropy at fixed circulation 
and energy [13]. In the viewpoint of Miller- Robert-Sommeria, this entropy describes the free 
merging of a system with two levels of vorticity ctq and a\ while in the other viewpoint, it 
describes the evolution of a forced system where the forcing has two intense peaks described 
by the prior x(o") = XoS{°~ ~ o"o) + XiK® ~ I 8 ]- 

The relaxation equations (12)-(13) can also be used as a numerical algorithm to construct 
stable stationary solutions of the 2D Euler equation. Indeed, Ellis et al. [7] have shown that 
the maximization of a functional of the form (10) at fixed energy and circulation determines 
a stationary solution of the 2D Euler equation of the form oj = f(ip), where / is monotonic, 
which is nonlinearly dynamically stable. Since the stationary solution of Eqs. (12)-(13) max- 
imizes S at fixed E and T (by construction), this steady solution of the relaxation equations 
is also a nonlinearly dynamically stable stationary solution of the 2D Euler equations (1). 
Thus, by changing the convex function C{oj) in Eq. (12), we can numerically construct a 
wide diversity of stable solutions of the 2D Euler equations. This is a potentially interesting 
procedure because it is usually difficult to solve the differential equation — Aif; = f(ip) directly 
and be sure that the solution is (nonlinearly) dynamically stable. These nonlinearly stable 
steady states can be an alternative to the statistical equilibrium state in case of incomplete 
relaxation, when the system has not mixed efficiently (non-ergodicity) so that the statistical 
prediction fails. In case of incomplete relaxation we cannot predict the equilibrium state but 
we can try to reproduce it a posteriori. 
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Finally, we have proposed in [11] to develop a phenomenological/effective statistical theory 
of 2D turbulence to deal with complex situations. The idea is that some types of entropic 
functional S[tJ] (in tJ-space) may be more appropriate than others to describe a given physical 
situation. For example, the enstrophy functional turns out to be relevant in certain oceanic 
situations [15] and the Fermi-Dirac type entropy in jovian flows [16, 9]. Certainly, other 
functionals of the same "class" would work as well for these systems. In addition, other 
classes of functionals S\p] may be relevant in other circumstances. Therefore, as a simple 
and practical procedure to describe a given system, we propose to pick a functional S[uJ] 
in the "class of equivalence" appropriate to that system and use it in the parameterization 
(12)-(13). We can thus describe the time evolution of the system on the coarse-grained scale. 
This approach is not completely predictive because we need to know in advance which type 
of entropy S[ZJ] describes best such and such situation. In practice, it must be determined 
by trying and errors (e.g. by comparing with oceanic data). But once a specific entropy 
has been found for a physical situation, we can work with it for different initial conditions 
specified by the robust constraints E and T (the effect of the Casimirs is reported in the 
chosen form of entropy SfuJ]). The idea is that the entropy S remains the same while E and T 
are changed. The problem is rich and non-trivial even if S has been fixed because bifurcations 
can occur depending on the control parameters E, T. This heuristic approach can be viewed 
as a simple attempt to account for the influence of the Casimirs while leaving the problem 
tractable. We use the fact that the Casimirs lead to non-standard (i.e. non-Boltzmannian) 
uj — ip relationships at equilibrium which are associated with non-standard forms of entropy 
S\uJ] in iU-space. We propose to fix the S'-functional depending on the situation. We do not 
try to predict its form, but rather to adjust it to the situation contemplated. This is based 
on the belief that some functionals S[uJ] are more relevant than others for a given system. 
Whether this is the case or not remains to be established. All the ideas presented here can 
be generalized to the case of quasi-geostrophic or shallow- water equations [9]. 
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